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Abstract 



By considering the Debye screening and the damping rate of glu- 
ons, the viscosity comcient of the quark-gluon plasma was evaluated 
via the real time finite temperature QCD in the relaxiation time ap- 
proximation at finite temperature and chemical potential. The results 
show that both the damping rate and quark chemical potential cause 
considerable enhancements to the viscosity coefficient of the hot dense 
quark-gluon plasma. 

PCAC number(s): 12.38.Mh, 12.38.Bx 

1 Introduction 

The calculation of the viscosity coeffcient in a hot and dense system is of 
interest both in the fields of high energy heavy ion collisions and astrophysics 
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[1-4]. The dissipative processes in a quark-gluon plasma(QGP) supposed to 
be formed in ultrarelativstic heavy collision can be described by the viscosity 
coefficient. 

Because of the estimated very short collision time of the ions in the heavy 
ion collisions, the equilibrium statistical approaches are probably inappropri- 
ate for an adequate description of such a process, and the dissipative phenom- 
ena are important, at least for the equilibration and the expansion phases 
of the system. In principle,the dissipative processes in heavy ion collisions 
should be described by non-equilibrium dynamical theory that is based on 
QCD. 

There are two methods to calculate the shear viscosity in QGP. The first 
one is by use of the relativstic kinetic theory. Starting from the Boltzmann 
equation, the viscosity coefficient r\ can be derived containing the transport 
cross section [2, 3]. The second one is from the kubo formulas[2]. For example, 
in the relaxation time approximation, the shear viscosity coefficient was found 
to be i] = cT 3 /a 2 s ln(l/« s ), the constant c was estimated to be 1.02, 0.28, or 
o.57,whereas a variational calculation gave c = 1.16 with two quark flavors 
in high temperature limit [2,3]. And Ref.3 also gave the result beyond the 
leading logarithm approximation, 

_ T 3 0.11 0.37 

On the other hand , basing on the kubo formula, Ref.2 gave r] < 2.6T 3 /a s 
and the lattice calculations obtained the value near the phase transition as 
< r]/T 3 < 9.5 [2]. 

It is known that the ensential nature of particles in a plasma at finite 
temperature is that they no longer have a infinite lifetimes, but have finite 
widths becaue of the collective effects. That is quite different from the be- 



havior at zero temperature [6,7]. In QGP Debye screening and damping of 
collective modes occur because of the interactions among particles and the 
heat bath. Therefore it is reasonable to take into account both the Debye 
screening and finite width effects when study the collective behaviors of the 
QGP, , which might also influence the transport processes of QGP. 

On the other hand, so far most all of the results of the shear viscosity 
cofncients in QGP are calculated at finite temperature without taking into 
account of the baryon density. Certainly by lattice simulation it is difficult to 
include the baryon density ,which is related to the baryo-chemical potential 
/it, through jj, = fib/ 3 and measures the deviation from the balance of quarks 
and and antiquarks. However at RHIC energies,there might still exist a 
considerable amount of stopping to /i/T ~ 1 — 2 [5] ,/i is the quark chemical 
potential. 

In this paper we will investigate the damping rate effect of the transverse 
gluon on the shear viscosity coefficient of QGP at both finite temperature and 
chemical potential quark chemical potential. We will show explicitly both 
the finite width and chemical potential cause considerable enhancements to 
the shear viscoisity. 

The paper will be organized as follows. In Sec. II we will get the shear 
viscosity of quarks and gluons via real time QCD by considering both the 
Debye screening and finite width effects at finite chemical potential. In Sec. Ill 
we will consider the running behavior of the coupling constant ane present 
numerical results. Finally we will summarize our results in Sec. IV. 

2 Calculation of the viscosity coefficent 

We will calculate the shear viscosity coefficient by using the relativistic ki- 
netic theory for a massless QGP in the so-called relaxation time approxima- 



tion[2,3]: 
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Vi = TzEiXi, (1) 

lo 

where is the energy density of particles of type i in the QGP and 
Aj is the mean free path, which is the inverse interaction rate: Aj = 
Considering the large angle scattering dominate the transport processes, we 
should multiply the interaction rate by a factor sin 2 0, 

T trsp . = J sin 2 0T t , (2) 

which is the so-called transport interaction rate, and is the scattering angle 
in the center of mass system. 

There are two equivalent ways of computing the interaction rates either 
using the matrix elements or from imaginary parts of the quark or gluon 
self-energies [3]. 

r, = -T^/mTr(^£)| P0=B . (3) 

Let us consider first the quark self-energy shown in fig.l, where we have 
included the screening effects and the damping rate effects by using the 
effective gluon propagator [8]. For hard quark momentum (< p q >~ T) 
it is sufficient to use the bare quark propagator and bare vertices. We will 
calculate the imaginary part of the self-energy at finite temperature and 
chemical potential using the Thermo Field Dynamics. For the hard particles 
one can show that the main contributions to the interaction rate comes from 
the soft momentum transfer range, i.e.g ~ gT in the weak coupling limit [3]. 

Using the usual Feynman rules of QCD to evaluate the diagram of Fig.l, 
one can obtain: 

Tr{^) = -9C f j— —— D = A-e, (4) 



where A^ u is the effective gluon propagator, in covariant gauge[9], 

A^q) = A^(q)A T (q) + B^A L (q) + D^A^q), (5) 

where 

A I \ X IX 9«?iu D (X 9m?0n q 2 , x QvOOs 

A^{q) = S^iPij - -Q2)0j V , Bp, = (V - -^-)q^{^0 ~ -^-) 

D».= qj f, \ = t- 2 , A T , L = ^— , (6) 
q 2 q 2 q 2 - H T ,L 

where £ is the gauge parameter. 

In the static limit the contribution of the hard thermal loops for hot dense 
QCD are[5], 



n L (g - 0, Q) = n 00 (g - 0, Q) = ml = g 2 T 2 (l + -± + JL £ %) (7) 



N f 1 ^ H 2 
v, v; — "(-o — y j- \l -r x 

for the longitudinal part and 

n T ( go ^o,Q) = o (8) 

for the transverse part of gluon polarization tensor. Substituting eqs.(7),(8) 
into eq(6),one gets 

' 1 ■ 1 (9) 



q z — rriQ q z 
In a thermal system, because of the thermal effects and the particles in- 
teractions, the particles (quasi-particles) will no longer have infinite lifetime, 
namely there are damping for the particles in QGP, which is an important 
feature of QGP. If one take into account the damping rate of the transverse 
gluon v T ~ Na s T[3,8], the transverse gluon propagator can be written as 
[6.7] 

By using eqs(5)-(10),one can cast eq(4) into: 
Tr( m = -iC, J ^ J£L ( AM,) - A M - (1 -f)^). (11) 



N 2 -l 



where C f - 2N . 

For finite temperature QCD, we will use the hot dense propagators for 
quarks and gluons in 2 x 2 matrix in TFD [9] . Among them the hot propagator 
of quark has the form 

iA n (p) = -tA 22 (p) = 2n(e(po)e^\ f (x p )+e(-p )e^/\(-x p ))5(p 2 )) 



iA 12 (p) = -iA 21 (p) = -27T ^e^(^(po)e^ /2 n / (^)+^(-Po)e^ /2 n / (-x p ))5(p 2 ) 

(12) 

where 

n f( x p) = , 1 , -, > x p = /5(po + 
e p + 1 

And for the gluons 



A n (g) = -iA£ 2 (g) = 5 - 2t™% 2 - m 2 )n B (g) = A^ + Aj, 

g z — TOq — «e 

Ai 2 (g) = 27ri<J(g 2 - m 2 )^^-^/ 2 , 
A n(?) = ^ e 7 — 2 - 2mlm- r^— ^(g), 

(g - ivt) 1 - q 2 (go - ivt) 2 - q 2 

An (go -> 0, g) = ^ ~ _ 2 7ri5(g 2 - z/ 2 )n B (g), 

A^ 2 (g) = 2/ m A T n B ( g )e-"l*l/ 2 
' ' ' 1 n B (g)e-^°l/ 2 , (13) 



«l?o| (go-^r) 2 -<5 2 (go + ^T) 2 -Q 2 

where is the Bose distribution function and for soft mementum g 

1 T 

nB ^ = e /%,| _ i ~ j^f ' ^ + ?) ~ n ^ p) - 



It is known that one can get the imginary part of the retard Green func- 
tion of the quark self-energy from the imaginary part of (1 — 2) type Green 
function, which is very convenient to calculate in TFD[10] 

e Pn/2 

Imt = —— JmEia, (14) 

sin 20 p+M 

where 

9( Po )e-^ + 9(-p )e^ 4 



sin < 



QXp/2 _|_ Q—Xp/2 



cos0 + _ 0(Po)e x ^ + 6(- Po )e-^ 4 

Then one can obtain 

r H D ~ 1 n 1 1 ens 2 f) 

ImTr(^)=4g*C f T / ^|b| 2 5((p+g) 2 )(^^-^^-(l-0^) 

J {2ii) z q l + v l q l + m 2 q l 

(16) 

Substituting eq(16),into eq(2), and assuming the transfer momentum q is 
soft: < g ~ gT < T, g/P ~ g we obtain 

- @ /* 2 d - o^JX^ - ^ - (i - 00(15,) 

which shows that the gauge dependent part is suppressed by the square of the 
coupling constant g 2 by evaluating the integration over q in < q ~ gT < T, 
we obtain: 

S 2 C,T r 2 , T 2 + m 2 ^, „ 2 v 

r <™=!^( mln ^^ + T ln ?^ ) ' (17) 

The transport interaction rate of gluon can be obtained in a similar way 
from the gluon self-energy, and it was found that 

N 

^trsp,g = -£T^trsp,q- (18) 



One can evaluate the quark and gluon energy density in QGP via the 
vaccum graphs in hot QCD both at finite temperature and chemical poten- 
tial: [1] 

By making use of eqs(l),(2),(19), and assuming p ~ 3T, v T ~ Na s T[8] 
we can get the shear viscosity coefficient of quarks for two flavors and gluons 
respectively 

21(l2ll T 3 H — 5tk M , 1 2rn _,_ (-, _ 2a s \\ 

h o2(4/3 + lna s + ln( ] 1 „ + a s ) + A{v)) ' 

16 / 7T 2 \rpZ ( 1 15a s \ 



„ _ 57LL3Q/ V 4tt / /21) 

19 a^/3+^)(-\na s + H^^ )+ a s )+A(u)y 

A{v) = —An N2c<X \ 5 ln(a 2 ^% A (22) 

In order to compare with previous calculations,one may come back to 
vanishing quark chemical potential case and not take the gluon damping rate 
into account, namely substituting A(v) = in eqs(20)-(22),then one obtains 
the shear viscosity coefficient of QGP with two quark flavors at vanishing 
quark chemical potential without including the finite width effect: 

T 3 ,0.25(1 - fg*) 0.73(1 - 15as ^ 



V ~ + w0.6+a s 4 ^ 1' ( 23 ) 



which agrees with the previous result in ref.3 beyond the leading logarithm 
approximation in the weak coupling limit except a factor 2 due to the defi- 
nition. At leading order and leading logarithm approximation, one can cast 
eq(23) into 

" = a98 ^7TK' (24) 



The coefficient c = 0.98 is remarkably close to the one found by Thoma 
c = 1.02 [3] and the one found by Baym et al.[2], the reason of the difference 
is we used < p q >~< p g >~ 3T other than < p q >~ 3.2T, < p g >~ 2.7T. 

3 Numerical analysis and discussion 

Before our numerical analysis, we consider the running behavior of the cou- 
pling constant g. In principle the running coupleing constant is governed 
by the renormalarization equations of QCD at finite temperture and chem- 
ical potential , which is still under investigation . Hereby we only use the 
following result [12]: 

a(T ' ^ = (33-2^)^2^^' (25) 
where is the QCD papameter, and we choose A s = O.lGeV as usual. And rif 
is the number of flavor, we will only consider two flavors. 

Substituting the running coupling constant eq(25) into egs(20)-(22), we 
obtain the experessions of the shear viscosity coefficients as functions of tem- 
perature and chemical potential. So we can plot the curves of the quark and 
gluon viscosity coefficients against temperature in Fig. 2(a) and Fig. 3(a) 
for different quark chemical potentials /j, — 0, 0.2, OAGeV with and without 
including the gluon damping effect in QGP respectively. From Fig.2(a) one 
can find the quark chemical potential effect may enhance the quark viscos- 
ity coefficient tremendously because the increase of quark energy density in 
QGP, while the chemical potential effect on the gluon viscosity is almost 
negligible. From Fig. 3(a) we find that the three curves of gluon viscos- 
ity coefficient against T at different chemical potential almost meet together 
both with and without taking into acount the gluon finite width. On the 



other hand the gluon damping effect causes large enhancements to the shear 
viscosities both for quark and gluon. The reason is that the damping effect 
has decreased their transport interaction rates, and hence increase their mean 
free paths. This qualitatively agrees with the analyses in ref.4 for the nuclear 
matter. 

Fig. 2(b) and Fig. 3(b) show the temperature dependence of the scaled vis- 
cosity coefficients rj/T 3 for different quark chemical potentials \i = 0, 0.2, OAGeV 
with and without including the gluon damping effect . Our result agrees with 
the Kubo formula result in ref.2 r//T 3 > 2.Q/a s at vanish baryon density, but 
our result is much bigger than the lattice estimate near the critical temper- 
ature in ref.2 . Fig. 2(b) and Fig. 3(b) showes at vanish chemical potential 
the curves of scaled viscoisity coefficients are no longer flat lines because of 
the running coupling constant. On the other hand, the tendencies are quite 
different between the scaled viscoisity of quark and gluon. The reason is 
the quark chemical potential plays different roles in them. Neverthless, the 
finite width still cause considerable enhancement to both the scaled viscosity 
coefficients 

Now we study how the shear viscosity coefficients vary with the coupling 
constant. We plot the curves of the scaled viscosity coefficients i]/T 3 against 
the coupling constant at a fixed temperature T = 300MeV and different 
quark chemical potentials /i — 0, 0.2, OAGeV with and without including the 
gluon damping effect in Fig. 2(c) and Fig. 3(c). The curves are almost above 
the Navier-Stockes Regime and in agreement with Thoma's result. Also 
the tendencies are quite different between the scaled viscoisity of quark and 
gluon via the coupling constant because of quark chemical potential plays 
different roles in them. However , the finite width still cause considerable 
enhancement to both the scaled viscosity coefficients. One interesting thing 



is that, although our calculations are valid for weak coupling limit g < 1, 
Fig. 2(c) and Fig. 3(c) show that the results can be extrapolated to higher 
coupling constant at fixed temperatur T = 300MeV, as predicted in ref.3 at 
least for quantities of energetic particles with E » T. 

4 Summary and Conclusions 

We have examined the quark and gluon shear viscosity behaviors by real 
time QCD both at finite temperature and baryon density in the quark gluon 
plasma in which both Debye screening and damping or finite width are re- 
garded as important underlying features. Also the running feature of the 
coupling constant is accounted for. We have seen that both the quark and 
gluon shear viscosity coefficients are increased substantiallydue to the ef- 
fects of the finite width of transverse gluons. Meanwhile all the scaled shear 
viscosity coefficients are increased in comparison with the previors results, in 
which the finite widtheffect was not accounted for. 

We have also examined the quark chemical potential effects to the shear 
viscosity coefficients in a straight forward way. It also enhances the shear 
viscosity coefficients of QGP consiberably, but the enhancements are smaller 
than that caused by the finite width effect. On the other hand, we have 
seen the tendencies of the quark viscosity coefficients, especially the scaled 
ones, differ quite substantially from that of the gluon. This is because the 
baryon density plays different roles in quark and gluon shear viscosity coef- 
ficients. 

We note that If we remove the quark chemical potential and finite width 
the results will be in good agreement with previous calculations [2,3,13]. 

We conclude that both the gluon damping rate and quark chemical po- 
tential effects have a considerable influence on the dissipative processes of 



QGP. This dissipation cannot be neglected in hydrodynamic descriptions of 
the expansion phase of QGP in ultrarelativistic heavy ion collisions. 
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Figure captions 

Fig.l The resummed one-loop quark self-energy graph 

Fig. 2 (a) The quark viscosity coefficient as a function of the temperture for 

/j, — 0, 0.2, 0.4 (bottom to top) without (thin solid lines in lower part ) 
and with (thick solid lines in upper part) including the gluon damping 
effect . 

Fig. 2(b) The scaled quark viscosity coefficient rj/T 3 as a function of the 
temperture for \i = 0, 0.2, 0.4 (bottom to top) without (thin solid lines 
in lower part ) and with (thick solid lines in upper part) including the 
gluon damping effect . 

Fig. 2(c) The scaled quark viscosity coefficient rj/T 3 as a function of cou- 
pling constant a at fixed temperature T = 300MeV for /i — 0, 0.2, 0.4 
(bottom to top) without (thin solid lines in lower part ) and with (thick 
solid lines in upper part) including the gluon damping effect . 

Fig. 3(a) The gluon viscosity coefficient as a function of the temperture for 
jj, = 0, 0.2, 0.4 (meet together) without (thin solid lines in lower part ) 
and with (thick solid lines in upper part) including the gluon damping 
effect . 

Fig. 3(b) The scaled gluon viscosity coefficient rj/T 3 as a function of the 
temperture for /i = 0, 0.2, 0.4 (bottom to top) without (thin solid lines 
in lower part ) and with (thick solid lines in upper part) including the 
gluon damping effect . 

Fig. 3(c) The scaled gluon viscosity coefficient rj/T 3 as a function of cou- 
pling constant a at fixed temperature T = 300MeV for /i = 0, 0.2, 0.4 



(bottom to top) without (thin solid lines in lower part ) and with (thick 
solid lines in upper part) including the gluon damping effect . 
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